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Predicting the progress of diffusively limited chemical reactions 
in the presence of chaotic advection 
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The effects of chaotic advection and diffusion on fast chemical reactions in two-dimensional fluid 
flows are investigated using experimentally measured stretching fields and fluorescent monitoring 
of the local concentration. Flow symmetry, Reynolds number, and mean path length affect the 
spatial distribution and time dependence of the reaction product. A single parameter A A, where A 
is the mean Lyapunov exponent and N is the number of mixing cycles, can be used to predict the 
time- dependent total product for flows having different dynamical features. 
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Chemical reactions in solution are enormously en- 
hanced by mixing, which brings the components into in- 
timate contact along lines or surfaces. While turbulent 
flows are effective in producing mixing, even laminar ve- 
locity fields can produce complex distributions of mate- 
rial and promote reaction through " chaotic advection" , 
in which nearby fluid elements separate exponentially in 
time. For fast reactions, the overall rate is determined by 
diffusion of reactants into the reaction zones, a process 
that is augmented by the stretching of fluid elements. 

The interplay of stretching, diffusion, and reaction has 
been investigated in various one-dimensional (ID) [l|, 0, 
13] and 2D models 0, d, d, 0, H in both open and closed 
domains. For irreversible, fast chemical reactions, nu- 
merical studies have attempted to relate product con- 
centration growth to the stretching properties of the flow. 
At early times, and assuming uniform stretching leading 
to exponentially growing reactant interfaces, the product 
growth is also expected to be exponential [i,[i3,[l!l. Sev- 
eral numerical investigations have shown that inhomoge- 
neous stretching should strongly affect local and overall 
reaction progress 0, d, [l2|- At later times, and using 
the passive scalar approximation in the presence of diffu- 
sion (valid for an infinitely fast reaction) , several different 
asymptotic functional forms for the product growth have 
been predicted, 0, S E, m • Boundary effects may lead 
to different regimes as time evolves 

However, there has been little opportunity to test the- 
oretical models because of the difficulty of measuring the 
stretching properties of experimental flows. Here, we use 
a recently developed method to measure stretching fields 
of two-dimensional time-periodic flows [ist . which quan- 
titatively describe the local finite time Lyapunov expo- 
nent field, and apply it to a fast reaction in order to ad- 
dress a fundamental question: How does chaotic mixing 
influence the spatial distribution of the reaction prod- 
uct, and the time dependence of product growth? While 
locally the reaction progress is related to the stretching 
properties of the flow, we show that the spatial average of 
the product concentration, which is proportional to the 
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FIG. 1: Streamlines computed from measured velocity fields 
at Re— 56 and path length p=2.5 for (a) disordered and (b) 
ordered magnet array configurations. The ordered case shows 
clear lines of mirror symmetry. Corresponding stretching 
fields for (c) disordered and (d) ordered fiows, showing the 
magnitude of stretching over an interval At= 1 period. Re- 
gions of strong stretching are localized in both cases. 



total quantity of product, depends on time and on the 
mean Lyapunov exponent in a way that is independent 
of the spatial structure of the flow (c.f. Fig. 4). We 
give a simple function describing this dependence. The 
experimental configuration, described in more detail else- 
where [161, is a thin fluid layer (1 mm thick) containing 
the reactants, that overlies a somewhat deeper conduct- 
ing fluid layer (3 mm thick). Chaotic flows are created 
by magneto-hydrodynamic forcing of the conducting fluid 
layer. A time-periodic electric current (with frequency 
/ in the range 10-140 mHz) passing through the lower 
layer, in the presence of an array of magnets, drives a 
time-dependent vortex flow that may be either spatially 
ordered or disordered depending on the magnet arrange- 
ment. The fluid is a 20% glycerol-water mixture (vis- 
cosity = 1.74 cP, density = 1.1 g/cm'^), about 10 cm x 
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10 cm in size, but only the central 8x8 cm is ima ged . 
The flow in the upper layer is nearly two-dimensional |17| 
and typical RMS velocities { U) are 0.05 to 0.7 cm/s. The 
Reynolds number, Re= pLU //u, based on the mean mag- 
net spacing (L=2 cm) is in the range 5 to 75. The path 
length parameter, p= U / Lf, which describes the mean 
displacement of a typical fluid element in one forcing pe- 
riod, is in the range 0.5 to 3.5. 

We investigate reactive mixing using an aqueous acid- 
base reaction NaOH + HCl ^ NaCl + H2O, or B^0+ 
+ 0H~ — > 2H2O, in the upper fluid layer. We write it 
schematically as A + B ^ 2P. The reaction is fast and 
second order. Its rate constant is fc=l.l x 10* s~^, 
and the reaction speed is characterized by the Damkoler 
number Da= kCoL'^/D, based on the ratio of the diffusion 
time scale {L? / D) to the reaction time scale (fcCb). Here, 
the diffusion constant D of either acid or base with its 
counter-ion is about 10^^ cm^/s, and Co=2.2 xlO^^ M 
is the initial reactant concentration. We note that Da 
is large (> 10^), so that the reaction is hmited by the 
diffusion fluxes toward the lines of contact. However, 
diffusion is enhanced by the stretching of interfaces. The 
relative importance of stretching and diffusion is given 
by the (Lagrangian) Peclet number Pe—L^\/D, where A 
is the mean Lyapunov exponent of the flow. Here, Pe is 
typically large and in the range 3.48 x 10^ to 1.02 xlO*". 

To determine fluid stretching, we first obtain high res- 
olution velocity fields using particle tracking methods 
described elsewhere 15|. Streamlines of the measured 



velocity fields for both ordered and disordered magnet 
array configurations at Re=56 (and p=2.5) are shown in 
Figs. l(a,b), respectively. The disordered flow has no 
spatial symmetry while the ordered flow has reflection 
and discrete translation symmetry along the coordinate 
axes. Regular (non-mixing) regions are found in the or- 
dered flow at a given Re, and in the disordered flow at 
low Re. 

Next, we use the velocity fields to construct displace- 
ment maps over a selected time interval At, and then 
compute stretching fields from the maps by differentia- 
tion [15]. This method, originally developed in numerical 
studies, provides high resolution stretching fields [H, [l^ . 
Examples are shown in Fig. l(c,d). The local stretching 
(S) measures the deformation of an infinitesimal circular 
fluid element located initially at {x,y) over the interval 
At. The local finite time Lyapunov exponent is defined 
as A=(log S)l At. Stretching fields computed over 1 pe- 
riod, corresponding to the disordered and ordered array 
at i?e=56 (p=2.5), are shown in Figs. l(c,d), respec- 
tively. Both fields show a wide distribution of stretch- 
ing values [2^ stronger at some locations than at others 
by a factor of 1000. In the ordered case (Fig Id), it is 
clear that the stretching is highly inhomogeneous, being 
much larger along lines passing through stagnation (hy- 
perbolic) points of the flow than in other regions. The 
strongest overall stretching occurs for the disordered ar- 




FIG. 2: (Color online) Normalized product concentration 
fields (P) vs. number of periods (TV) at p=2.5. (a) -Re=37 
(disordered), 7V=10; (b) Same as (a) for A^^SO; (c) _Re=56 
(disordered), iV=10; (d) Same as (c) for Af=30; (e) 7?e=56 
(ordered), A^^IO; (f) Same as (e) for A''=30. Fully reacted 
regions appear in red, as shown in the color bar. Unreacted 
regions are shown in black. 



ray (Fig. Ic), due in part to the lack of spatial symmetry 
of the velocity field. We also find stronger stretching for 
flows with larger p at a given Re and magnet array con- 
figuration. The mean Lyapunov exponent A is computed 
by averaging the local values over the image domain. 

We now show in Fig. 2 the reaction of initially seg- 
regated aqueous solutions containing reactant A (acid) 
and reactant B (base) at i?e=37 and i?e=56 for the dis- 
ordered flow, and also at iZe=56 for the ordered flow. 
Initially, a solid barrier (the dotted line in Fig. 2a) sepa- 
rates the acid, which contains a pH sensitive fluorescent 
dye, from the base. The barrier is lifted and the reaction 
is observed over varying numbers N of cycles. (A^^IO 
and A^=30 are shown in the upper and lower rows of Fig. 
2.) 

We define normalized concentration fields for the 
acid A, base B, and product P, as follows: A—A/Aq, 
B—B/Bf), and P~P / Pfinai, where ^0 and Bq are the 
initial reactant concentrations, and P final is the product 
concentration in the fully reacted state. Using separate 
calibration experiments, we determine A{x, y, t) from the 
local light intensity. From the conservation of material 
expressed in the statement A -I- B ^ 2P, it can be shown 
that averaged over the entire cell {P)~1-2{A) , where 
{A)=l/nJ2^=i and n is the number of pixels in the 
image. Though this relation is not accurate locally due 
to advection and diffusion, we use P= 1-2^4 as an approx- 
imation to display the local normalized product concen- 
tration field. As time evolves, the interplay of stretching, 
diffusion, and reaction creates a complex pattern, with 
regions of high (red) and low (dark) normalized product 
concentration. The development of convoluted interfaces 
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FIG. 3: (Color online) (a) Spatial averaged normalized prod- 
uct concentration (P) as a function of time (or A^) and Re 
for several different flows (D=disordered, 0=ordered, all at 
p=2.5). (b) (P) vs. mean Lyapunov exponent ( A) for various 
N (disordered flows for various Re) . Lines are added to guide 
the eye. Product growth is slow for small but accelerates for 
large at later times. 



is evident early. At longer times, systematic differences 
between the different flows become clear. For i?e=37 
(disordered array), the reacted regions are spatially ex- 
tended, but large unreacted regions remain even after 30 
periods. This flow departs less from time- reversibility, 
a necessary condition for mixing [lij . due to its lower 
inertia than does the i?e=56 case. Even in the reacted 
regions, the product concentration is small at lower Re. 

Product formation is also affected by the extent of reg- 
ular (non-mixing) islands, which are favored by the mir- 
ror symmetry of the ordered flow (Fig. 2e,f), but that 
also occur for the disordered flow at the lower Re (Fig. 
2a,b). These are regions of low stretching where the in- 
terface between reactants grows approximately linearly 
in time, and the local reaction rate is less than in re- 
gions of high stretching. Such isolated regions are not 
visible at in Fig. 2(c,d) (higher iJe=56, disordered flow), 
where the accumulated product is higher than it is in the 
ordered flow for the same Re and reaction time. These 
qualitative differences show that the product distribution 
is substantially affected by the flow pattern. 



Next we explore the evolution of the spatial average 
of the normalized product (P) (or the total quantity of 
product) as a function of time for flows having differ- 
ent dynamical features. Fig. 3(a) shows that the high 
Re and p disordered flow produces the fastest product 
growth, while lower i?e, p or spatial order substantially 
reduces the growth rate. The reason is that irreversibil- 
ity due to inertia (large Re), large particle displacement 
(p), and the absence of barriers to transport (disorder) 
are required for rapid mixing. For all flows, the initial 
product growth is roughly linear rather than exponen- 
tial in time. This is especially clear for low Re and or- 
dered flows, where there are substantial low stretching 
regions. Although the simplest theories suggest expo- 
nential early growth llj this is not apparent in the data, 
most likely due to the wide distribution of local stretch- 
ing rates first pointed out in Ref . [3] ■ A numerical study 
in 2D fiows Q also finds that non-chaotic regions lead to 
sub-exponential initial product growth. 

The average normalized product (P) is most usefully 
parameterized by the mean Lyapunov exponent A rather 
than Re or p. In Fig. 3(b), we show the variation of (P) 
as a function of A at different numbers of periods N for 
disordered flows with different Re {p=2.5). The growth 
with A is approximately linear for small A but accelerates 
for large A at the later times. 

Although the concentration patterns are complex, the 
evolution of (P) (at least after an initial transient) turns 
out to be a function only of A'^ and A. In Fig. 4, we plot 
(P) vs. AA^ for the disordered flows at various Re and 
p. We find that product concentration curves collapse 
onto a single master curve with no adjustable parameters. 
This is a surprising result because even though the flows 
have the same magnet array configuration, they possess 
different degrees of time reversibility and mean particle 
displacements per cycle, and regular regions occur to a 
different extent in the various flows. In the insert to Fig. 
4, we show that this scaling behavior can be extended to 
ordered flows as well. Remarkably, the rescaled ordered 
flow data for different Re fall onto the same master curve 
as the disordered array cases, which include various Re 
and p. 

The dependence of product concentration on and time 
has been recently considered theoretically 2l| using a 
simplified model for the reaction interface. The authors 
suggest that after a short time (i<l/A) the normalized 
concentration should be described by 



(P) = 1 - exp{-aXN) 



(1) 



where a is a constant. A similar relationship is also found 
in a numericaj_ study using periodic boundary conditions 
for 2D flows [71]. This result, shown by the dotted line in 
Fig. 4 for the best fitting value of a describes our data 
adequately for AA^ < 300. However, for larger values, 
(P) grows more rapidly than this form predicts. 
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FIG. 4: (Color online) Dependence of (P) on the product of 
the mean Lyapunov exponent (A) and period TV for disordered 
flows at various Re and p. Insert: Early behavior for both 
ordered and disordered flows (at various Re and p). The 
dotted and solid lines correspond to Eq.l (a=0. 00263) and 
the empirical fit described in the text, respectively. 



for uniform stretching. Spatial symmetry, particle dis- 
placement, and the departure from time reversibility af- 
fect the reaction progress (and A). However, the normal- 
ized product concentration up to the fully reacted state 
approximately follows a single master curve with XN as 
the independent variable, where A is the measured mean 
Lyapunov exponent and N is the number of mixing cy- 
cles. This result allows quantitative prediction of the 
total product as a function of time for several chaoti- 
cally mixing flows with different dynamical features such 
as spatial symmetry and the extent of departures from 
reversibility. This scaling is not expected to apply to tur- 
bulent flows, but may be appropriate for a variety of fast 
reactions that are stirred by chaotic advection. 

We thank Z. Neufeld and G. Haller for fruitful discus- 
sions, and T. Tel and G. Karolyi for communicating their 
theoretical results to us. T. Shinbrot provided helpful 
comments on the manuscript. This work was supported 
by NSF DMR-0405187. 



Other theoretical works show a double exponential 
growth, or a more complex spatio-temporal behavior de- 
pending on the reaction rate [13] or the presence of 
boundaries Finally, one might consider a solution 

of Floquet form consisting of a sum of several exponen- 
tials 22], since the flow is periodically driven. However, 
none of these forms is found to fit the data over the entire 
range of AA^. 

We can fit the entire data set (ordered, disordered, and 
different Re and p) to a single empirical function of the 
form (P) = 1 - exp[-aXN - l3(XNf], where a=0.0018 
± 1.9 X 10"'' and /3= 6.2 x 10"^ ± 2.5 x lO"'^ are con- 
stants. (However, the constants and could depend on the 
initial shape of the reaction interface, a variable we did 
not explore.) Note that this equation (solid line in Fig. 
4) fits the data over the entire range reasonably well. For 
the lower values of AA^, the quadratic dependence is weak 
and we recover Eq.(l). The more rapid approach to sat- 
uration at late times might be due to the transport of 
initially segregated reactants to regions where they can 
react. A similar transport effect was demonstrated by 
Voth et al. fl^ for passive mixing. Note that we cannot 
reach the fully reacted state for A < 3. 

In conclusion, we study the interplay of stretching, 
diffusion, and fast reaction experimentally, with direct 
measurement of the stretching fields. The strong spatial 
heterogeneity of stretching, with a distribution spanning 
many decades, [3 [2^ causes the early product growth to 
deviate strongly from the exponential behavior expected 
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